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Abstract, We employ Magnetohydiodynamic (MHD) simu- 
lations to examine the creation and evolution of plasma vor- 
tices within the Earth's magnetosphere for steady solar wind 
plasma conditions. Very few vortices form during intervals 
of such solar wind conditions. Those that do remain in fixed 
positions for long periods (often hours) and exhibit rotation 
axes that point primarily in the x or y direction, parallel 
(or antiparallel) to the local magneto spheric magnetic field 
direction. Occasionally, the orientation of the axes rotates 
from the x direction to another direction. We compare our 
results with simulations previously done for unsteady solar 
wind conditions. By contrast, these vortices that form dur- 
ing intervals of varying solar wind conditions exhibit dura- 
tions ranging from seconds (in the case of those with axes 
in the x or y direction) to minutes (in the case of those 
with axes in the z direction) and conveet anti sunward. The 
local-time dependent sense of rotation seen in these previ- 
ously reported vortices suggests an interpretation in terms of 
the Kelvin-Helinholtz instability. For steady conditions, the 
biggest vortices developed on the dayside (about 6 Re in di- 
ameter), had their rotation axes aligned with the y direction 
and had the longest periods of duration. We attribute these 
vortices to the Hows set up by reconnection on the high- 
latitude magnetopause during intervals of northward Inter- 
planetary Magnetic Field (IMF) orientation. This is the first 
time that vortices due to high-latitude reconnection have been 
visualized. The model also successfully predicts the principal 
characteristics of previously reported plasma vortices within 
the magnetosphere, namely their dimension, flow velocities, 
and durations. 


Keywords. Magneto spheric physics (MHD waves and 
instabilities) 


1 Introduction 

The solar wind— magneto sphere interaction can be steady or 
unsteady. Unsteady processes, like transient vortices at the 
magnetospheric boundary, may dominate the interaction by 
transferring mass, energy and momentum into the magneto- 
sphere. Many processes have been suggested for the creation 
of such vortices, like pressure pulses, magnetic reconnec- 
tion and the Kelvin— Helmholtz (KH) instability (Sibeck et 
aL, 2002: Shi et ah, 2005: Mima, 1995: Kivelson and Chen, 
1995). The latter is the most commonly cited for the cre- 
ation of the vortices and occurs when surface waves are ex- 
cited on the outer magnetospheric boundary by the veloc- 
ity shear provided by magneto sheath plasma flowing past 
the magnetosphere. 

Observations from single spacecraft, in conjunction with 
theoretical models and simulations, have frequently been 
used to explore and understand the KH instability (Ogilvie 
and Fitzenreiter, 1989: Chen et ah, 1993: Otto and Fairfield. 

2000: Hasegawa et ah, 2006: Takagi et ah, 2006: Fairfield 
et ah, 2007; Collado- Vega et ah, 2007). These studies ex- 
plored the parameters favoring development of the instabil- 
ity. Some also addressed the transport of mass and momen- 
tum by the instability. Using multi- space craft observations. 
Ha sega wa et ah (2004) employed Cluster observations to 
provide direct evidence for rolled-up vortices on the flank 
magnetopause during an interval of generally northward 
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Interplanetary Magnetic Field (IMF). Hasegaw r a et aL (2005) 
compared these observations with the output from three- 
dimensional magnetohydro dynamic simulations for the solar 
wind's interaction with the magnetosphere. More recently, 
Palermo et al. (2010) presented KH vortex signatures de- 
tected by two of the five THEMIS spacecraft near the post- 
noon magnetopause. 

Although the KH instability may be the best known in- 
stability capable of creating vortices near the magnetopause, 
other processes may also generate vortices. For example, 
plasma flows into the reconnection region, is accelerated, and 
flows outward. During intervals of southward IMF orienta- 
tion, we expect reconnection on the day side equatorial mag- 
netopause. Magneto spheric plasma on closed magnetic field 
lines flows into the reconnection regions, while reconnected 
plasma on newly opened magnetic field lines flows antisun- 
ward along the magnetopause (Dungey, 1961). In the case of 
southward IMF near local noon, the flows would be in the 
meridional plane with an axis nearly in the y direction. 

By contrast, during intervals of strongly northward IMF, 
we expect IMF field lines to be appended to the magneto- 
sphere via high latitude reconnection poleward of both cusps. 
There is flow into the reconnection site within the magneto- 
sphere and the newly reconnected magnetic field lines move 
anti sun w-ard along the flanks of the magnetosphere (Song and 
Russell, 1992; Kessel et al., 1996; Song et al, 1999). This 
again results in flows with axes nearly in the y direction. Both 
kinds of reconnection result in shears and large-scale plasma 
convection vortices (Saunders, 1989). 

Indeed, vortices with solar wind properties very differ- 
ent to those predicted for the KH instability have definitely 
been observed, Zhang et al (2011) reported vortical flows 
near the Low Latitude Boundary Layer (LLBL) in conjunc- 
tion with flux transfer events (FTEs) during southward IMF, 
They discuss how anti sunward moving FTEs transfer mo- 
mentum to the LLBL plasma through a viscous interaction. 
This interaction drives antisunward flows in the LLBL along 
the sides of the FTEs, increasing the thermal pressure ahead 
of it, diverting the plasma normal to the magnetopause and 
then allowing it to return to locations where the pressure 
is low. forming a vortex. Flow vortices can also be related 
to field-aligned currents with possible corresponding auroral 
signatures (Sitar et al,, 1998). 

Col la do -Vega et al. (2007) examined the output of MHD 
simulations driven by real solar w ind input. They found that 
vortices generated under northward IMF were more preva- 
lent on the dawns ide than on the dusks ide, but were sub- 
stantially less ordered on the dawns ide than on the dusk- 
side. The scale sizes of most of the vortices w^ere large, up 
to 10 /?£. Rotation axes w^ere closely aligned with the z di- 
rection in solar magnetic (SM) coordinate system. The vor- 
tices rotated preferentially clockwise on the dawns ide, and 
counterclockwise on the dusk side. These characteristics sug- 
gested an interpretation in terms of the KH instability. 


A question arises as to whether such vortices w r ould de- 
velop during intervals of steady solar wind conditions. To 
isolate cause and effect, this paper examines the magne- 
to spheric response for fixed solar wind conditions, namely 
fixed density, temperature, and velocity. We compare and 
contrast the case for a nominal 360 kins 1 solar wind speed 
with one for a higher 700 kins"" 1 solar wind speed. In each 
case, w r e rotate the IMF direction from southward to north- 
ward as a step function 15 min into the run. We study how r 
file IMF direction affects the growth and polarization of the 
vortices, and how fast the vortices develop depending on the 
velocity shear applied to the boundary. Based on our previous 
study, we predict that the vortices develop near the equatorial 
plane, with axes in the z axis direction and having a sense of 
rotation that depends on location (clockwise on the dawn- 
side and counterclockwise on the dusk side). Some vortices 
did develop with these characteristics, but most of them had 
very different characteristics. 

Section 2 describes the simulation and methods used for 
this study Section 3 describes the vortices found near the 
magnetopause during fixed nominal solar wind conditions 
(case 1) and for high-speed solar wind conditions (case 2). 
Section 4 discusses and compares the vortices found in this 
study with those in the previously reported simulation for 
dynamic solar wind conditions (Collado- Vega et al.. 2007). 
Closing remarks constitute Sect. 5. 


2 Methodology 
2.1 MHD simulation 

For this study w r e initiated the MHD simulation with a 
step function in the IMF direction; 15 min being southward, 
— 5 nT, and then turning northward, +5nT, for tW'O hours. 
Two cases, one w r ith a nominal speed solar wind (360 km s -1 ) 
and another w r ith a high speed solar wind (700 kins -1 ), w r ere 
nm using the global MHD Block-Adaptive-Tree-Solarwind- 
Roe-Upwind- Scheme (BATS-R-US) code (Pow r ell et al., 
1999; Gombosi et al., 2004: Tdtli et al., 2012) developed at 
file Center for Space Environment Modeling at the University 
of Michigan. BATS-R-US is part of the Space Weather Mod- 
eling Framework, and w r as available to us by the Community 
Coordinated Modeling Center (CCMC). 

BATS-R-US global magnetosphere (GM) solves 3-D 
MHD equations in finite volume form using numerical meth- 
ods corresponding to Roe's approximate Riemann solver. It 
is coupled wdth the ionospheric electrodynamics solver that 
solves Poisson's equation. The code has a block- adaptive 
grid composed of three- dimensional rectangular blocks ar- 
ranged in varying degrees of spatial refinement levels. It 
uses a grid with maximal resolution of 0.25 in the in- 
ner magnetosphere, the magneto sheath, the magnetopause 
and plasma sheet. The magnetospherlc MHD part is attached 
to an Ionospheric potential solver that provides electric 
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Fig. 1. Example of a vortex which would not be obvious with the naked eye : but that is correctly found with our algorithm. The colors 
represent the velocity magnitude in log scale and the arrows represent the flow vectors. The left panel is 2-D in the z = 0 cut plane on the 
dusks ide where the vortex is located. The right panel is the 3-D view of the vortex rotating approximately on the x axis. 


potentials and conductances in the ionosphere from magne- 
tospheric field-aligned currents (Ridley et ah, 2004). It uses 
local time stepping, which allows each grid cell to converge 
toward steady state as fast as possible, limited only by the 
local numerical stability limit. 

The input parameters are transformed to geocentric solar 
magneto spheric (GSM) coordinates, and the Earth's mag- 
netic field is approximated by a dipole with updated axis 
orientation and co -rotating inner magnetospheric plasma or 
with a fixed orientation dining the entire simulation run. 
The orientation angle is updated according to the time simu- 
lated or a fixed axis position can be specified independently 
from the time interval that is simulated. In our simulation 
the dipole had a fixed orientation of 0 degrees for the en- 
tire simulation run, and a resolution of 0.4 Re was used on 
the day side, which is small compared to the typical size of 
Keivin-Helmholtz vortices (y 3—10 /?e). 

In a previous study (Collado- Vega et aL, 2007). statis- 
tics of vortices found near the ecliptic plane on the mag- 
netopause flanks were presented using the Lyon-Fedder- 
Mobbary (LFM) MHD simulation (Fedder et aL, 1995a, b; 
Lyon et aL, 2004). The simulation was performed in the solar 
magnetic (SM) coordinate system, allowing for the tilt of file 
Earth's magnetic dipole relative to the solar wind flow direc- 
tion to be included. Outflow^ conditions are imposed on the 
downstream boundary. Elsewhere, external boundary condi- 
tions were specified using ACE solar wind data propagated 
appropriately to the front and cylindrical sides of the grid. 


2.2 Vortex identification and visualization 

Although the concept of a vortex is familial - , a rigorous defi- 
nition has been elusive. Lugt (1995) suggested defining a vor- 
tex as “the rotating motion of a multitude of material particles 
around a common center”. When viewing a vector field of ve- 
locity, we can use tins definition to recognize voitices. How- 
ever, quantitative measures are necessary to progress beyond 
labor-intensive manual identification. A measure of persis- 
tency is also intuitively part of the criteria, hence the link with 
coherent structures in a flow, A mathematical quantity, vor- 
ticity, is defined as the curl of the velocity vector (a> = V x t?). 
It is a vector quantity whose direction is the axis of rotation, 
and is related to the local angular - rate of rotation. Vorticity 
magnitude and isosurfaces are used to detect vortical struc- 
ture successfully in some flows. Hie criterion of a local pres- 
sure minimum at a vortex arises from swirling flow arising 
from a balance between centrifugal force and pressure gradi- 
ent. but tills is not a sufficient condition and is problematic 
in 3-D. 

An alternative way of viewing the flow is by its stream- 
lines, the instantaneous curves that are everywhere tangent to 
the direction field defined by the velocity vectors. This leads 
to another qualitative definition for vortex identification: “A 
vortex exists when instantaneous streamlines mapped onto 
a plane normal to the vortex core exhibit a roughly circu- 
lar or spiral pattern, when viewed from a frame of refer- 
ence moving with the center of the vortex core” (Robin- 
son. 1991). There are some difficulties associated with im- 
plementing this criterion, in particular it is not Galilean in- 
variant. which means that whether or not the streamlines are 
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Fig. 2. Time duration of vortices for case 1. The x axis represents the number of vortices and the y axis represents the time in hours for each 
of the 4 quadrants. The blue X and Y arrows are just to show the separation in each of the four quadrants, X positive being in the sunward 
direction and Y positive in the dusk ward direction. Hie vortices rotating in the clockwise direction are represented in red color, and the ones 
rotating counterclockwise are represented in blue. 


closed depends very much on one's frame of reference. Jeong 
and Hussain (1995) regard this as a requirement, although 
others, ourselves included, have had success in implementing 
this criterion. 

Mathematical techniques adapted from topology have 
proven useful in elucidating structure in vector fields (Perry 
and Chong, 1994). Critical points are points in the flow field 
where the streamline slope is indeterminate and the velocity 
is zero relative to an appropriate observer. These points can 
be used to intuit important aspects of die vector field topol- 
ogy. These concepts may be more familiar in their use for 
magnetic field topology and visualization (Cai et al.. 2001). 
In pursuit of vortices, the velocity gradient tensor is of par- 
ticular filter est. Let Vv denote the velocity gradient tensor. 
Perry and Chong (1987) propose defining a vortex core as a 
region with complex eigenvalues of Vv. The premise is that 
complex eigenvalues indicate closed streamlines in a refer- 
ence frame moving with the point, V v can be written as a 
sum of a symmetric pail S = | (Vv + Vv T ) and an antisym- 
metric part £2 = “(Vu- Vv 1 ). Then Q = |(|£2 2 -\S 2 ), the 
second invariant of the velocity gradient, is a measure of the 
shear strain rate and vorticity. Hunt et al. (1988) used the 
criterion of positive Q to define eddies, along with low T pres- 
sure. As Hunt et al, (1988) observed, for most vortices tills 


criterion is equivalent to complex eigenvalues of the velocity 
gradient tensor. All these methods w r ork well for some fluid 
flow r s, but our experience thus far with simulation data for the 
Earth’s magnetosphere indicates that these measures are not 
sufficient to uniquely identify space plasma vortices. 

Our current algorithm takes a multi- staged approach to lo- 
cate vortices in 3-D flows. It requires only velocity compo- 
nents at each point in the domain. The second invariant of the 
velocity gradient tensor. Q, is used to identify points where 
the rotational strength exceeds the strain rate. The points that 
pass this test (2 > 0) are then subjected to further scrutiny. 
The vorticity vector is used to define the £-axis of a local 
coordinate system, and the velocity field is projected onto 
a plane perpendicular to this vector. Local streamlines are 
computed. A point is classified as a vortex if a small neigh- 
borhood exhibits both closed streamlines and a velocity mag- 
nitude minimum at the center. Tins algorithm is success- 
fill at finding the vortex center, but is not immune to false 
hits. More work is required to fine tune the search criteria 
and characterize both vortex size and strength. One method 
of screening for alleged eddies due to numerical effects is 
to duster the vortices into classes by spatial proximity. The 
class size can be used as a filter to eliminate isolated points 
that are more likely to be numerical artifacts. It is far more 
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Table 1 . Solar wind plasma input parameters. 


Parameters 

Case 1 

Case 2 

V x 

“360 kms -1 

— 700 kms -1 

Vy 

— 28 kms - 1 

-28 kms -1 

V z 

—38 kms" 1 

-38 bus" 1 

IMF B x 

2nT 

2nT 

IMF By 

1 nT 

1 nT 

IMF B z 

— 5/+5 nT 

-5/+5nT 

IMF B 

5.10 nT 

5.48 nT 

IMF clock angle 

168.7 deg 

168,7 deg 

Density 

7ncc~ 1 

3 nec” 1 

Temperature 

40 000 K 

40 000 K 


likely that a genuine vortex is part of a vortex tube that has a 
spatial extent. Figure 1 shows an example of a vortex which 
would not be obvious with the naked eye, but is correctly 
found with our algorithm. At first, the visualization tool only 
shows 2-D cut planes of a box with dimensions of x = —25, 
25 /?e, y = —20. 20 Re The vortex algorithm provides the 
location from which to begin so the locations of interest are 
identified faster, and can be visualized in the 3-D view. The 
left-hand image shows a cut plane from a data simulation 
(cropped, near the tail at the magnetosphere) with z — 0. The 
vectors show the x and y components of the velocity vec- 
tor and the colors represent the velocity magnitude in a log- 
arithm scale. The right-hand image shows the 3-D velocity 
field for the data within the white square in the left-hand 
image. 

Although our vortex algontlmi is very successful at iden- 
tifying individual vortices, it is less successful at grouping 
vortices both temporally arid spatially. In order to determine 
whether individual vortices were really part of the same vor- 
tex, for example, persisting at one location for an extended 
time or convecting tail ward over time, two additional meth- 
ods were employed. The MHD Explorer tool discussed in 
Collado- Vega et aL (2007), was used to visualize both indi- 
vidual vortices and groups of vortices. The human eye is ca- 
pable of discerning collections and patterns for which com- 
puter algorithms do not yet exist. Secondly. Excel spread- 
sheets were used to gr oup similar characteristics (axes of ro- 
tation, rotation direction, diameter, etc.), determine persis- 
tence hi time, and to plot voitex locations after each vor- 
tex was carefully studied by the naked eye. The resulting 
collected vortices are discussed in the following sections. 

3 Vortices under fixed solar wind conditions 

In Table 1 we describe two sets of input parameters for the 
MHD simulation. For each case we fixed the solar wind ve- 
locity (Vjc- V y , Vy, IMF (Z?*, By. B z ), IMF clock angle, 
density, and temperature. Each simulation ran for a total of 
2h and 15 mill. For both cases the IMF B z component was 



Fig. 3, Dawnside vortex rotating clockwise with its rotation axis 
aligned in the z direction (pattern 1 ) that developed during south- 
ward IMF within the first 15 min of the simulation run. 


initially southward for the first 15 min. and then turned north- 
ward with the same magnitude (5 nT) for the rest of the simu- 
lation. In case 1 the x component of the solar wind speed was 
nominal with a value of —360 kms -1 . The transverse compo- 
nents had significantly lower values (V^ = -28 kms -1 and 
V z — —38 kms -1 ), the density was 7ncc -1 and the temper- 
ature had a value of 40 000 K. The IMF transverse compo- 
nents had small values: B x = 2nT and B v = 1 nT. and die 
IMF clock angle was 168.7 degrees. For die lhgh solar wind 
speed case (V x = 700 kms L ), die density had a lower value 
of 3ncc 1 , in keeping with observations of typical solar 
wind trends. The x and y components of the IMF and the 
clock angle were die same as in case 1. 

Using the automated vortex detection algorithm, spread- 
sheets. and the MHD Explorer visualization tool with the 
output parameters of the BATS-R-US simulation runs, a total 
of 33 vortices were found for the two cases with fixed solar 
wind conditions; 18 vortices were found for the slower solar 
wind speed case, and 15 vortices were found for die faster 
solar wind speed case. The specifics are discussed hi the next 
two sections. 

3.1 Nominal solar wind speed: case 1 

During the 2 h and 15 min simulation run, vortices developed 
throughout the simulation in two distinct patterns: (1) vor- 
tices with their rotation axis in the z direction, and (2) vor- 
tices with their rotation axis in either the x or y direction. The 
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Fig. 4, Dusk me lit side vortex that shows how some vortices could be rotating oil the x axis and then twisting into another direction. 


first pattern includes vortices that developed during the first 
15 min with southward IMF, vortices that developed early on 
when the IMF turned northward, and vortices that developed 
much later in Hie simulation. For the first pattern, vortices 
developed along both the dawn and dusk flank; those on the 
dawnside rotated in the clockwise direction, while those on 
the duskside rotated preferentially in the counterclockwise 
direction (of four cases, three rotated in the counterclockwise 
direction). For the second pattern, vortices were not as well- 
ordered and could twist so that part of a vortex would have a 
rotation axis in the x direction, and another part of the same 
structure would have its rotation axis in the y or z direction. 
Their rotation directions u r ere both clockwise and counter- 
clockwise, not dependent on location. These formed more 
than tliiity minutes after the simulation began when condi- 
tions were steady. Figure 2 shows vortices foimd in each of 
the four quadrants, along with their duration in the simula- 
tions. The dusk magnetotail had many more vortices than 
the other quadrants. Vortices on the dayside ty pically lasted 
longer than on the night side. All of the vortices are described 
in more detail below'. 

During the first 15 mm of the simulation when the IMF 
w r as southward, only 2 vortices developed, as shown in 
Fig. 2, quadrants 2 and 3 (Q2, Q3). One vortex developed 
on the dawn tail flank, x — -0.79 Re, y — - 13 Re* arid the 
other developed on the dusk tail flank, x = —2.76 /? E , y = 
12.72 Re< Both of the vortices fit the first pattern, with their 


rotation axis closely aligned to the z direction. In Fig. 3 we 
show' the vortex that developed on the dawnside. Its rotation 
axis is closely aligned parallel to the positive z axis and it 
rotates in the clockwise direction, where the direction of ro- 
tation for all the vortices is always chosen when viewed from 
the positive coordinate axes. The vortex that developed at the 
duskside (not shown) rotated in the counterclockwise direc- 
tion. No vortices developed on the dayside during the 15 min 
of southward IMF (Q1 and Q4 in Fig. 2). 

Ten minutes after the IMF turned northward, one vortex 
fully developed (Q3 in Fig. 2), This was along the dusk tail 
flank, with a rotation axis in the z direction and counterclock- 
wise rotation. This vortex lasted 30 min and slowiy connected 
down tail. On the dawn side tail flank (Q2), a pattern 1 vortex 
developed and then dissipated, with another one developing 
at about the same position 10 min later and then converting 
down tail. This vortex also had its axis in the z direction, 
but rotated clockwise. After about one and a half hours of 
simulation time, a few r more vortices that fit the first pattern 
developed in the dusk tail. These were short lived and found 
in the vicinity of vortices of the second pattern. 

Pattern 2 vortices were found in all four quadrants. On the 
dayside they formed after about an hour of simulation time 
and persisted for tire rest of the simulation. The dayside pat- 
tern 2 vortices had their rotation axes in the y direction; there 
w^as one clockw ise and one counterclockwise in both dayside 
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Fig. 5. Inner dusk vortex rotating counterclockwise with its axis 
closely aligned with the x direction (pattern 2), 

quadrants (QI and Q4 in Fig. 2). These vortices formed in die 
z~ —5.6 Re cut plane. 

On the nightside, the pattern 2 vortices were more likely to 
be oriented in the x direction but could twist into another di- 
rection. An example is shown in Fig. 4. Here a dusk nightside 
vortex formed around one hour and 30 min of the simulation 
run with one part rotating oil the x axis, and then turning into 
the z axis. On the dawnside the two pattern 2 vortices rotated 
in the clockwise direction. On the dusk side pattern 2 vor- 
tices were found with both rotation directions, clockwise and 
counterclockwise. In Fig. 5 a nightside dusk vortex is rotat- 
ing counterclockwise with its axis closely aligned with the x 
direction. 

Table 2 lists the characteristics of the vortices that devel- 
oped for the two patterns in each quadrant. This table shows 
that the day side vortices were significantly larger than the 
nightside vortices f 6 Re on the day side compared to 2—\ Re 
on the nightside. The table summarizes the characteristics 
that w ere described above. 


3.2 High solar wind speed: case 2 

For the high solar wind speed case, case 2, a total of 15 
vortices formed. As for case 1, vortices also developed in 
the same two distinct patterns. The first pattern formed with 
the vortex axes in the z axis direction, and the second pat- 
tern formed by the vortex axes in the x-y direction. Fig- 
ure 6 shows the vortices that formed under case 2 in all four 
quadrants with their time duration. 

During the first 15 min of the simulation, when the IMF 
was completely southward, only two vortices formed. They 
both developed at the beginning of the simulation (Fig, 6 Q2 
and Q3), One developed at the dawn tail flank (Q2) and the 
other at the dusk tail flank (Q3), The dawnside vortex rotated 
clockwise, and the dusks ide vortex rotated counterclockwise. 
Both vortices had their rotation axis aligned to the z axis. As 
in case 1, no vortices developed oil the day side region for 
the southw ard IMF. After the 1 5 min of southward IMF had 
passed, only one more vortex developed with its axis closely 
aligned to the z axis (pattern 1). It developed at around thirty 
minutes before the end of the simulation at the dawn tail flank 
and moved towards the boundary (— z) very slowly until the 
end of the simulation. 

Pattern 2 vortices started developing after 35 min of the 
IMF being northw r ard. These vortices developed at the dawn 
and dusk dayside with their rotation axis closely aligned to 
the y axis (Fig. 6, Ql and Q4), They w^ere the bigger sized 
vortices developed for this case, with about 6 Re of diameter. 
They are also seen until the end of the simulation run, lasting 
for about 1 h and 25 min. Some vortices also developed at 
the dawn and dusk nightside region. However, more vortices 
fanned at the dawnside for case 2. These vortices had the 
rotation axis closely aligned to the x axis. They w r ere mostly 
about Re in diameter, and lasted around 30-45 min. No 
preference in rotation w as noted for the vortices of pattern 2. 
These vortices' durations can be seen in quadrants 2 and 3 
in Fig. 6. 

Also, two vortices in pattern 2 developed simultaneously 
at the same x-y location but with a different z coordinate, 
and they rotate in the opposite sense. In Fig. 7 the vortices 
ar e shown having approximately the same x-y location, but 
with different z coordinates on the dusk dayside. One of 
the vortices rotates clockwise, and the other rotates counter- 
clockwise. Vortices also appeared to have symmetry in w hich 
dawn and dusk vortices would develop simultaneously on the 
magnetopause flanks with different rotation directions. Fig- 
ure S presents a 2-D cut plane with z = — 5.6 Re, where a 
dawn and a dusk vortex developed simultaneously rotating 
counterclockwise (dusk) and clockwise (dawn). 

Table 3 lists the char acteristics of the vor tices fonned for 
case 2. It show r s that the dayside vortices w^ere bigger (6 Re) 
than the biggest ones fonned at the nightside (3-4 Re)- Also, 
it show r s that for this case the majority of the vortices were 
fonned at the dawnside, not at the deskside Like in case 1 . 
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Fig. 6. Same as Fig. 2, but for the vortices found for case 2 (higher solar wind speed). The vortices rotating clockwise are represented in red 
color and those rotating counterclockwise are represented in blue. 


Table % Case 1: four quadrant vortex statistics. 


Pattern 


Dawn dayside 



Dawn night side tail 

# 

Axis 

Rotation 

Size 

# 

Axis 

Rotation 

Size 

1 

0 




2 

Z 

cw 

3-1 R e 

2 

2 

7 

ccw.cw 

6 Re 

2 

X 

cw 

3—4 % 

Pattern 


Dusk dayside 



Dusk nightside tall 


# 

Axis 

Rotation 

Size 

U 

Axis 

Rotation 

Size 

1 

0 




4 

z 

ccw. cw 

2-4 « E 

2 

2 

7 

ccw,c\v 

<S*E 

6 

X 

ccw.cw 

2-4 R e 


4 Discussion and comparison 

This study explored vortices developed in an MHD simula- 
tion with fixed solar wind conditions, while a previous study 
(Collado- Vega et ah, 2007) looked at dynamic solar wind 
conditions shown in Fig. 9. Using the new search algorithm 
described in Sect. 2 we re-examined data from this previous 
study, but only times for which the conditions w ere similar to 
the current study, with a dynamic pressure of about 5— 10 nPa. 
The magnetic field, however, was greater. 10-15 nT for one 
interval and 20-25 11 T for the second interv al (cases C and E 
from Collado -Vega et al., 2007). The most significant differ- 


ences betw een the vortices formed with dynamic solar wind 
conditions and the ones formed under fixed solar wind condi- 
tions are rotation (preference and axis), convection, location, 
time of duration and other characteristics as described below-. 

4.1 Rotation 

In the current study, pattern 1 vortices for both case 1 and 
case 2 had axes aligned with the z direction and rotated pref- 
erentially clockwise on the dawnside. and counterclockwise 
on the duskside. For both cases, a total of 4 pattern 1 vortices 
w-ere found on the dawn mghtside. from which 3 rotated in 
the clockwise direction, and a total of 5 pattern 1 vortices 
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Fig. 7, Vortices on the dusk day side. The vortices formed almost at the same x-y location, but with a different z coordinate. Tins is in the 
y = 7 cut plane. 


Table 3. Case 2: four quadrant vortex statistics. 


Pattern 


Dawn dayside 



Dawn nights i de tail 

# 

Axis 

Rotation 

Size 

u 

Axis 

Rotation 

Size 

1 

0 




2 

Z 
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3 % 

2 

2 

y 

ccw,cw 
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6 

X 
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3-4 % 

Pattern 
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# 

Axis 

Rotation 

Size 

U 

Axis 

Rotation 

Size 

1 

0 




1 

z 

ccw 

3 He 

2 

2 

y 

cw,ccw 

fi* E 

2 

X 

cw 
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were found on the dusk nights ide. from which 4 rotated in 
the counterclockwise direction. Pattern 2 vortices had rota- 
tion axes mostly aligned with the x or y axis. Oil the day- 
side, case 1 and case 2 pattern 2 vortices had rotation axes 
aligned with the y axis, but on the liightside their rotation 
axes were aligned with the x axis. No preference in rotation 
direction w r as evident for pattern 2 vortices. In the previous 
dynamic study, vortices developed with rotation axes in the 
z axis. In this study these w r ould be classified as pattern 1 
vortices. However, in the re-examination of the two similar 
cases with our new automated algorithm* 71 vortices w : ere 
found with rotation axes aligned with the x or y direction. 


These would be classified as pattern 2 vortices in our current 
study. Figure 10 shows a dawn vortex of 30 March 2002 at 
00:40 UT corresponding to the higher magnetic field case of 
the dynamic conditions for the previous study. Tliis vortex 
has its rotation axis closely aligned to the x axis. 

In the previous study, pattern 1 vortices rotated preferen- 
tially clockwise on the dawns ide and counterclockwise on 
the dusks ide as w r as seen for the pattern 1 vortices in our cur- 
rent study. This behavior u r as attributed to the polarization 
of surface waves that could be driven by KH instability at 
the magnetopause (Collado- Vega et aL* 2007). In our current 
study, more pattern 2 vortices developed than pattern 1 for 
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Fig. 8, Dawn and dusk vortices found simultaneously rotating counterclockwise (dawn) and clockwise(dusk). These vortices came in pairs 
with the ones shown in Fig, 7. The rotation direction is determined by following the arrowheads from the streamlines from the positive y axis. 
The vortices are located in the same z = “5.6 /?£. but the dawn dayside on Jt = 3.8 y = -7.2 R e. and the dusk dayside on x = 2.6 R e, 
y = 6.8 R 


both case 1 and case 2. The axes of these vortices tended to 
lie parallel (or antiparallel) to the magneto spheric magnetic 
field, and they also developed during times of stable solar 
wind conditions contrary to our previous study where vor- 
tices formed under dynamic solar wind conditions. 

Figure 11 shows the rotation polarity for all the vortices 
that exhibited rotation axes either parallel (counterclockwise 
when viewed from the dominant magnetic field axis) or an- 
tiparallel direction (clockwise when viewed from the dom- 
inant magnetic field axis) to the dominant magnetic field 
component. The sense of rotation therefore follows the right- 
hand rule. Dayside vortices follow a pattern with the ones on 
the duskside rotating parallel to the local magnetic field, and 
those on the dawnside rotating antiparallel. We understand 
tliis rotation in tenns of the global convection pattern and re- 
connection; plasma on closed magneto spheric magnetic field 
lines moves sunward through the dayside magneto sphere, 
while plasma on open magnetospheric magnetic field lines 
moves antisunward. No definite pattern was seen for the vor- 
tices that fonn on the nights ide. However, more vortices ro- 
tate parallel to the magnetic field oil the dawnside than on 
the duskside. More statistics need to be considered to reach 
a final decision about the nightside vortices. The black curve 
shows the magnetopause location as calculated by the BATS- 


R-US model for the last time step in case 2. This shows that 
all the vortices were found inside the magnetosphere. 

4.2 Convection 

In the current study there were several pattern l vortices (z 
direction axis) that convected tail ward. These vortices oc- 
curred mostly during times when the solar wind conditions 
had just changed: at the beginning of the simulation and 
soon after the IMF switched from southward to northward. 
In case I there were four convecting vortices: two each for 
southw ard and northward IMF, on each tail flank. For case 2 
only tw T o convecting pattern 1 vortices formed, one on each 
dawn and dusk flank during southward IMF, In the previous 
study pattern 1 vortices were consistently swept antisunward 
by convection, becoming larger in diameter and changing ve- 
locity as they traveled along the flanks (Collado- Vega et ah, 
2007), As an example, in Fig. 12, a clear progression in vor- 
tex development is shown. The vortex is traveling antisun- 
ward along the dusk flank. Within only four minutes (02:46— 
02:49 UT). the vortex core travels approximately 12 Re down 
the flank. In the re-examination of the tw r o cases from the pre- 
vious study, the pattern 2 vortices that were found appeared 
to be convecting down tail: however, their time duration w T as 
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2002 

Fig, 9, Average values for die high-speed solar wind stream from 2S March to 5 April 2002. The shaded region shows the nine hour time 
S ame of the vortices searched hi Collado-Vega et al. (2007). 


ill terms of seconds. The longest duration vortices found only 
lasted for about 2 min. These characteristics will be described 
m more detail below. 

Most of pattern 1 vortices convected tailward for our cur- 
rent study. We found a few others that did move towards the 
boundary in the negative z direction. Once steady state con- 
ditions were established, convection disappeared. We found 
three more vortices for case 2 pattern 2 that moved but hi the 
negative z direction also. The vortices that convected tail- 
ward did it more slowly than vortices in the previous study, 
taking up to 30 min to convect down tail out of the simu- 
lation box. The rest of the vortices formed for case 1 and 


case 2 did not convect. They had fixed positions, and lasted 
longer periods of time than those that did convect. These re- 
sults suggest that it takes a change in solar wind conditions to 
generate vortices that convect tailward, but the change can be 
a sharp step function. The pattern 2 vortices found in the re- 
examination appeared to convect tail ward, but because most 
of them lasted just seconds, and since we are looking at one 
minute plasma date files with our visualization tool, it is dif- 
ficult to be completely sure about in which direction they are 
moving. 
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Fig, 10, Dawn vortex of the previous study case from Collado- Vega 
et al. (2007) of 30 March 2002, 00:40 UT. that was found using our 
automated search algorithm with its rotation axis closely aligned 
with the x axis. 


4,3 Location 

In our current study, pattern I vortices tended to be near the 
equatorial dawn and dusk flanks on the night side, for both 
connecting and non-convective vortices - This is also true for 
the vortices found in our prev ious study. No pattern 1 vortices 
developed on the dayside, Pattern 2 vortices developed in all 
4 quadrants; however, more developed on the dusk night side 
for case 1 and more on the dawn nightside for case 2. Pat- 
tern 1 vortices generated under northward IMF for the pre- 
vious study (Collado- Vega et aL 2007) were more preva- 
lent on the dawnside than on the duskside, and were sub- 
stantially less ordered on the dawnside than on the duskside. 
This is also true for the pattern 2 vortices found in the re- 
examination cases. By contrast, Taylor et al. (2012) found 
vortices with a dawn-dusk asymmetry, with 71 % on the dusk 
flank suggesting preferential growth in this region. 

There may be ionospheric signatures cone spending to the 
magneto spheric vortices. We would expect to observe the 
ionospheric signatures at the footprints of the magnetic field 
lines that pass through the center of the magneto spheric vor- 
tices. We mapped an elliptical-shaped set of points in one 
of the longest lasting day side vortices found in case 1 (and 
shown in Fig. 13) to a small area in the early afternoon (^ 
12:00-1 8:00 LT) northern and southern high-latitude (75°) 
ionosphere, as indicated by the white crosses in Fig. 14, The 
solid black line in the figure marks the boundary of the open 


magnetic field line region. The ellipse in the magnetosphere 
maps to a small vortex at the convection reversal boundary 
in the Northern Hemisphere as indicated by the flow arrows 
at the circle marked by the white crosses. No corresponding 
vortex is seen in the Southern Hemisphere. 

4.4 Duration 

For our current study vortices for both cases had long pe- 
riods of duration. However, pattern 2 vortices had substan- 
tially longer periods of duration than pattern 1 vortices, par- 
ticularly on the dayside. These pattern 2 vortices lasted for 
hours (as shown in Figs. 2 and 6). Pattern 1 vortices for the 
previous study (not shown) lasted minutes before they con- 
vened tailward and started to dissipate. Those found in the 
re-examination (pattern 2) lasted just seconds with, some of 
them surviving up to a few minutes. This is shown in Figs. 15 
and 16. The vortices are separated in the same four quadrants 
as we did for case 1 and case 2 of the fixed cases, but now 
for the two cases of magnetic field magnitude. Most of the 
vortices only lasted only a minute until they completely dis- 
appeared. Those that lasted more than a few r minutes were 
located in the nightside region. They are not as stabilized as 
those formed with rotation axes on the z direction for the dy- 
namic cases, which lasted several minutes before starting to 
dissipate. 

In the current study, most vortices did not convect, instead 
they rotate with axes parallel (or antiparallel) to the magnetic 
field. This could be the cause of slower convection vortices 
with longer periods of duration. For case 1 of the 18 total 
cases, 12 are rotating in the x or y direction. Of those 12,11 
are located off the equator at z > ±2.5. For case 2, w'e have a 
total of 15 cases; 12 of them are rotating in the x or y direc- 
tion and 10 of those 12 are located at z > ±2.5, In Fig. 17 we 
plotted the magnetic field lines near the rotation axis of the 
vortex shown in Fig. 5 using the Space Weather Explorer tool 
from the Community Coordinated Modeling Center at NASA 
Goddard Space Flight Center. The magnetic field lines are 
shown in red and the plasma flow' is shown with yellow lines. 
The vortex rotation axis is parallel to the magnetic field lines. 
The color background represents the plasma density. 

4.5 Other characteristics 

Other characteristics for the fixed conditions cases vortices 
are that they appear to have symmetry, and they would de- 
velop in pahs at approximately the same x-y location, but 
w r ith a different z coordinate. Some vortices would develop 
simultaneously on the daw'll and dusk side for the same time 
frame, This is shown in Fig. 8 where tw'o vortices with al- 
most the same size, rotation velocity and rotation axis de- 
veloped at the dawn and dusk side. How r ever, the dawnside 
vortex rotates counterclockwise, and the duskside vortex ro- 
tates clockwise. Figure 7 shows an example of tw r o duskside 
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Fig. 11. X-Y plot of the vortex rotation directions referred from the right-hand rule following the local magnetic field dominant component. 
The ones on the dayside followed a pattern; the ones on the duskside rotated parallel to the field while die ones on the dawn side rotated 
antiparallel. The black curve represents the magnetopause location as calculated by die BATS-R-US model for the last time step in case 2, 


vortices that formed almost with the same x-y coordinates, 
but with totally different z coordinates. 

For the dynamic solar wind conditions cases (previous 
study), 47 % of the total development of vortices formed un- 
der low speed solar wind conditions, before the onset of the 
high speed stream (see Collado- Vega, et ah, 2007, for refer- 
ence). This is also the case with the static solar wind con- 
ditions (current study), where 1 8 vortices formed under low 
solar wind speed (360 km s and 15 vortices formed dur- 
ing high solar wind speed (700 km s -1 ). However, it is ev- 
ident in both studies that it is necessary to have dynamic, 
stronger, faster and driven solar wind conditions to create 
faster tailward-convecting vortices. 

We attribute dayside vortices like the one shown in Fig. IS 
to high- latitude reconnection during northward IMF. The si- 
multaneous appearance of such vortices in both the Northern 
and Southern Hemispheres (see Figs. S and 7) suggest simul- 
taneous reconnection is occurring poleward of both cusps. 
As suggested by for example, Dungey (1963) and also il- 
lustrated by Reiff (1984), magneto spheric plasma flows sun- 


ward toward the high- latitude reconnection site, reconnects, 
and then flows antisunward and poleward along the magne- 
topause. In steady state, the sunward flows fonn the equator- 
ward edge of the vortices, while the antisunward flows fonn 
the poleward edge. Flows are very weak at the center of the 
vortices, whose locations remain fixed for the steady- state 
solar wind conditions in these rims. 

The magnetosheath/boundary layer flow drags the closed 
field lines that have recently been added to the magneto- 
sphere anti sunward to the magnetotail, where they slip az- 
imuth ally into the magnetotail and ultimately contract earth- 
ward within the center of the magnetotail. Vortices in the 
magnetotail with axes in the x direction result from the shear 
between newly appended magnetic field lines folding into the 
magnetotail at higher latitudes and the same magnetic field 
tines collapsing earth 1 ward deeper within the magneto tail at 
lower latitudes. 

In previous studies, waves and vortices were also ob- 
served on the same scale size and time duration as the ones 
foimd in this study. Couzens et al. (1985) observed waves of 
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Fig, 12* Voit ex progression of 30 March from 02:46 to 02:49 UT from Co 11a do- Vega et al. (2007). Within only 4 min the vortex core travels 
approximately 12 /?£ down the flank. The velocity of rotation increases significantly along the flank, from 20 to 60 kins - * . while its traveling 
velocity decreases when the vortex becomes fully ordered around 02:48 UT, This is not typical of the vortices found in the present study. 

Table 4. Previous observations. 


Authors 

Satellite 

Size 

Period 

Hones etaL(L97S, 1981) 

ISEE 1 and 2 

20-40 « E 

Up to 9 h 

Sckope et al. (1981) 

ISEE 1 and 2 

1 «E 

Few minutes 

Couzens et al (1985) 

ISEE 1 and 2 

11000 km 

Several hours 

Chen et al. (1993) 

ISEE 1 and 2. IMP S 

Wavelength 15 /?£ 

5 min 

Otto and Fairfield (2000) 

Geotail 

Approximately 2 fig 

2.5 min 
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Fig. 13. Cut plane for y = 6.80 for a case 1 dusk day side vortex. The color represents the velocity magnitude and the arrows the flow 
direction. The white marks represent the selected data points were field lines were selected to map to the ionosphere. 
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crosses mark the location wher e die vortex maps into the ionosphere at about 75 degr ees latitude between 12:00 and 1 S:00LT. A vortex is 
actually present in the Northern Hemispher e ionosphere, but not in the Southern Hemisphere ionosphere. 


approximately 1 1 000 km that lasted for several hours using 
ISEE I and 2 data, Chen et ai (1993) observed waves w ith 
an average period of about 5 mm with a wavelength of the or- 
der of 15 % using data from ISEE 1, 2 and IMP 8. Hones et 


al. (1978, 1981) published the first evidence for large vortices 
with a size of 20 R e near the plasma sheet that had periods of 
hours, Borovsky et al. (1997) showed v x -v y autocorrelations 
in the plasma sheet that indicate vortices with periods in the 
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Fig. 15. Time duration of vortices for case C (dynamic cases) with a 10-15 nT magnetic field magnitude. The vortices last only a few minutes 
before starting to dissipate. This figure only shows the duration of the vortices and not their rotation dir ection. 


range of 6-1 S min. Our events dimensions fluctuate from 2 
to 6 /?e, and in both cases (case 1 and 2) we have vortices 
that have long periods of duration, some in terms of hours. A 
summary of previous observations similar to those found in 
tliis study are presented in Table 4. 

Foullon et al. (2008) stated that the source region of KH 
waves varies in position due to the influence of the global 
magnetosphere geometry and that not much attention was be- 
ing applied to this factor In tliis study w r e analyzed the influ- 
ence of two different solar wind speeds and a step function of 
the magnetic field on the magnetosphere, and we do suggest 
that the vortices' development depends on their local posi- 
tion and magnetosphere geometry. Vortices with fixed posi- 
tions were found that could be the result of having steady 
solar wind conditions. 

To be certain that the results shown in this study are not 
model dependent, a steady solar wind condition mil with the 
same solar w r ind parameters used in case 1 above was made 
using the Lyon-Fedder-Mobbary global MHB code, Tliis 
w r as the code used in the study of Collado- Vega et al. (2007). 
We used our automated approach and the MHD Explorer vi- 
sualization tool to find vortices. Vortices exhibited the same 
characteristics (eg,, rotation axis and rotation direction) at 
similar locations to those described in this study, which used 
the BATS-R-US MHD code. Although the LFM code gen- 
erates more vortices and these vortices convect somewhat 


faster than those in the BATS-R-US runs, the vortices de- 
scribed in this study were found using both codes, meaning 
that they are not an artifact of ail individual code. 


5 Conclusions 

In this study, fixed solar wind conditions were simulated to 
examine the creation and evolution of vortices in a global 
MHD model for the Earth’s magnetosphere. A total of 18 
vortices were found for the case with nominal solar wind 
speed, 360 kins and 15 vortices for the cases of higher 
solar wind speed, 700kms _l . By contrast. Collado-Vega et 
al. (2007) presented statistics fox a total of 304 vortices found 
near the ecliptic plane on the magnetopause flanks, using 
simulated MHD data driven by real solar wind conditions, 
Foullon et al. (2010) stated that the Ke 1 vin-Helmholtz in- 
stability w r as more likely to occur for wave propagation in 
the direction of high flow shear. We do not see more vor- 
tices for intervals of higher solar wind velocity, which may 
result in a higher shear oil the flanks, and most of them do 
not move antisunward with their rotation axes in the z direc- 
tion. Therefore, we conclude that most of the vortices found 
in the current study w r ere not caused by the KH instability. 
Consequently, they could be associated with transients in the 
solar wind-magnetosphere interaction. The largest dayside 
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Fig. 16. Time duration of vortices for case E (dynamic cases) with a 20-25 nT magnetic field magnitude. In tins case the vortices last between 
2 and 5 min before starting to dissipate. The figure only shows their time duration. 



Fig. 17. A 3-D overview of the vortex in Fig. 5 using the Space Weather Explorer tool. The vortex rotation axis is parallel to the magnetic 
field lines. The background color represents the density, the yellow lines represent the flow and the red lines are the niagnetospheric field 
lines. 


www.anii-geophys.iiet/31/1463/2013/ 


Ann. Geophys., 31 * 1463-1483, 2013 



1480 


Y. M, Collado- Vega et ah: Vortices obstipations and comparison 



Fig. 18, Vortex developed at the dawn day side with the rotation axis closely aligned to the y axis (left). The same vortex can be encountered 
after 1 h and 16 min of simulation (right). This vortex has nearly a fixed position for more than an hour, but it actually moves closer to the 
magnetopause boundary in the sunward direction. 


vortices seen in this study were the result of high-latitude 
reconnection happening under northward IMF. This is the 
first time that vortices formed under high- latitude recon- 
nection have been visualized. Some of the vortices seen in 
the nights ide region were also formed by the reconnected 
magneto spheric field lines that moved along the magneto- 
sphere sides and got pushed back to the magnetosphere in 
the far tail region. The vortices were formed by the shear 
between the transition of the anti sunw ard flows of recently 
reconnected field lines and the sunward flow s of magneto- 
spheric closed field lines. 

The results of a detailed intercomparison of the new' results 
in this paper for steady solar wind conditions and the results 
previously obtained for time- varying solar wind conditions 
reveal: 

- Pattern 1 (z axis) vortices for both fixed and time- 
varying solar w r ind conditions rotate preferentially 
clockwise on the dawns ide and counterclockwise on 
the duskside. This pattern is consistent with expecta- 
tions for vortices driven by surface weaves on the mag- 
netopause. perhaps from the KH instability. By contrast, 
the sense of rotation seen in pattern 2 ( x and y axis) vor- 
tices does not depend on location, 

- Pattern 1 (z axis) vortices that form during intervals of 
time- varying solar wind conditions convect antisunward 
in the manner expected for surface waves. The majority 
of pattern 1 (z axis) vortices that form during intervals 
of steady solar wind conditions also convect antisun- 
ward, albeit more slowly. By contrast, pattern 2 ( x or y 
axis) vortices seen during interv als of steady solar wind 


conditions, like those shown in Figs. IS and 19, remain 
hi place, sometimes dissipating. Figure 19 shows a dusk 
flank vortex that remains fixed for about 30 min of the 
simulation before starting to dissipate. The Fig. 18 vor- 
tex. w ith an axis in the y dir ection, is part of the global 
convection pattern. 

— Vortices developed near the equatorial dawn and dusk 
flanks on the nightside for both cases in the present 
study. However, more vortices developed at the dusk 
nightside for case 1 (low r solar wind velocity) and more 
vortices developed at the dawn nightside for case 2 
(high solar wind velocity). The reasons for these differ- 
ences remain unclear. Pattern 1 and pattern 2 vortices 
generated under northward IMF for the previous study 
W'ere more prevalent on the dawnside than on the dusk- 
side. and were substantially less ordered on the dawai- 
side than on the duskside. 

— Vortices formed under fixed conditions had substan- 
tially longer periods of duration than vortices formed 
under dynamic conditions. Vortices formed under fixed 
conditions had periods of duration in terms of hours. 
Those formed under dynamic conditions formed and 
dissipated in terms of several minutes. This difference 
can be related to the vortices convection velocity for 
each case. For our current study, vortices rotate either 
parallel or antiparallel to the magnetic field, 

— As for the dynamic cases where most of the vortices 
formed before the onset of the high speed stream, more 
vortices for fixed conditions formed under the nominal 
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Fig. 19. This dusk flank vortex lasts for more than 30 min. All four boxes are exactly ten minutes apart but have the same spatial dimensions. 
The vortex core is not swept tailward. it remains fixed and it starts to dissipate around 02:12 UT (last box). 


case of sola] wind speed. For case 1 we found a total 
of IS vortices* but for case 2 a total of 15 vortices were 
found. 

- Seventy- one vortices were found rotating on the x or y 
direction for the dynamic cases (previous study) using 
our new search algorithm. However, these vortices have 
very short durations compared to the ones described in 
our previous study, and the two cases presented in tills 
paper, 

- We followed field lines of an elliptical-shaped set of 
starting points for one of the clays ide vortices found 
in case 1, and found a vortex in the Northern Hemi- 
sphere ionosphere around 12:00-18:00 local time at 75 
degrees latitude. These vortices form on closed mag- 
netic field lines and rotate either parallel or antiparallel 
to the field, which suggests that their footprints couple 
to the ionosphere. 


- The results shown in this paper are not model depen- 
dent. 
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